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ABSTRACT 

Recombination of the primordial helium plasma (Hell^Hel, z ~ 1500 — 3000) is 
considered. This process has an effect on the CMBR anisotropy and CMBR spectrum 
distortion. In this work an influence of neutral hydrogen on kinetics of Hell— >HeI 
recombination is investigated in the frame of the standard cosmological model. It is 
shown that small amount of neutral hydrogen (10 -5 — 10~ 2 of total number of hydrogen 
ions and atoms) leads to acceleration of Hell— >HeI recombination at z < 2000 and at 
z < 1600 quasi-equilibrium Hell— >HeI recombination (according to the Saha formula) 
becomes valid. 
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1 INTRODUCTION 



The Hell— >HeI primordial recombination has an ef- 
fect on formation of the cosmic microwave background ra- 
diation (CMBR) anisotropy (e.g. Hu et al., 1995; Seager et 
al., 1999) and CMBR spectrum distortion (Lyubarsky and 
Sunyaev, 1983; Fahr and Loch, 1991; Dubrovich and Stol- 
yarov, 1997; Wong et al., 2006). The relative difference of 
CMBR anisotropy power spectrum corresponding to Saha 
Hell— +HeI recombination and corresponding to Hell— >HeI 
recombination by Seager et al. (1999) has a level of up to 
5% for multipoles I ~ 1500-3000 (Seager et al., 2000). Such 
changes can be measured in future experiments (e.g. Planck 
and others). 

Kinetics of Hell— >HeI recombination has been consid- 
ered in a number of papers (Boschan and Biltzinger, 1998; 
Novosyadlyj, 2006; Wong and Scott, 2006; and references 
therein). The most used model of helium recombination is 
the model suggested by Matsuda et al. (1969) and devel- 
oped by Seager et al. (1999). A significant modification of 
this model was suggested by Dubrovich and Grachev (2005) 
who took into account recombination through Hel ortho- 
states. 

In the model by Seager et al. (1999) the effect of neu- 
tral hydrogen on Hell— >HeI recombination was neglected 
although in a number of previous papers (e.g. Hu et al., 
1995; Boschan and Biltzinger, 1998; and others) authors 
pointed out that this effect may be essential. This contra- 
diction made us investigate the effect of neutral hydrogen 
on Hell— >HeI recombination again and derive Hell— >HeI re- 
combination kinetic equation which takes this effect into ac- 
count. 



2 ADDITION TO THE RECOMBINATION 
MODEL 

The main addition to the Seager's et al. (1999) model is 
taking into account an interaction of resonant Hel quanta 
with neutral hydrogen. This interaction is described by the 
following elementary processes: 

1. Hel resonant photon ionizes hydrogen atom with emis- 
sion of electron which has a kinetic energy more than 6 eV: 
H + 7 -> H+ + e" . 

2. The emitted electron loses the kinetic energy during 
electron-electron collisions with thermal (E £t th = 3/2fcsT < 
0.7 eV) electrons of the primordial plasma. The typical time 
of this process can be estimated by the following formula 
(e.g. Spitzer, 1978): 



4TvN e e 4 ln{Am e vi/3k B T) ' 



A = 



2e 3 



(k B T) 3 



(1) 



where m e is the electron mass, v e is the electron velocity, 
N e is the free electron concentration, T is the temperature 
of medium. Under considered conditions the value t ee is less 
than 10 3 s. 

3. The hydrogen ion recombines with a thermal electron and 
emits a photon with energy ~ (13.6 + ksT) ~ 14 eV: H + + 
e~ — > H + 7. The typical time of this process is larger than 
10 9 s. 

Comparison of the typical times of electron-electron col- 
lisions and hydrogen recombination shows that between the 
ionization and the following recombination of hydrogen the 
energy of the non-thermal electrons is divided among the 
thermal electrons of the primordial plasma. 
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In the chain of elementary processes mentioned above 
Hel resonant photons disappear and do not keep the popu- 
lations of Hel excited states. 

The rate of hydrogen ionization by Hel resonant 2p — > 
Is photons Jh [cm~ 3 s _1 ] determines the rate of disappear- 
ance of these photons. It was shown by Seager et al. (2000) 
that this rate is negligible in comparison with the total rate 
of 2 1 p *-* l 1 s transitions (i.e. the difference of rates of for- 
ward 2 x p — > l 1 s and backward l 1 s — > 2 1 p transitions) of 
Hel due to escape of resonant photons from line profile be- 
cause of cosmological redshift Jnei,bg (subscript b denotes 
2 1 p state of Hel, subscript g denotes l 1 s state of Hel, fig. 
1) q In this work we have found that the rate Jh was un- 
derestimated by Seager et al. (2000). To show this fact we 
considered the ratio of the rates Jh and Jnei,b g - 

The rate Jh is given by the expression 



Jh — caH{v bg )N 1 , bg NHi 



(2) 



where c is the speed of light, an^Vbg) is the hydrogen ion- 
ization cross section at frequency v>b g of the 2 1 p — > Us tran- 
sition of Hel, Nhi is the concentration of hydrogen atoms 
in ground state, N Jt h g is Hel 2 1 p — » Is resonant photon 
concentration given by the formula 



N- 



Av D 87T^, 



2k B T 



where r\b g is the number of photons per mode 2 1 p — > Us, 
Aud is the thermal width of line 2 1 p — > l 1 ^, rune is the 
mass of helium atom. 

The rate JHei, bg [cm~ 3 s~ 1 ] is given by the equation 

JHeI,b g — Al g NHeI,b (4) 

where Al g is the effective coefficient [s^ 1 ] of 2 1 p <-> l 1 s tran- 
sitions due to photon escape from the line profile because of 
cosmological redshift (superscript r denotes "redshifting" ) , 
NHei.b is the concentration of Hel atoms in state 2 1 p. The 
effective coefficient of transitions is given by the formula: 



At 



bg 



c 3 g b N H ei, g 



(5) 



where H = #o\/^a + Sl m (l + z) 3 + ^rei(l + z) 4 is the 
Hubble constant as function of redshift (parameters Ho, Oa, 
Qm and Q re i described in tab. 1), gj, = 3 is the statistical 
weight of 2 p state, Nnei, g is the concentration of Hel atoms 
in I s state. 

Optical depth of absorption in the line Hel 2 1 p — > l 1 s 
is much larger than unity, therefore r/ bg is satisfied the fol- 
lowing expression 

NHei.b ,„s 

Vbg ^ JTT (6) 

g b ^HeI, g 

Using equations © - ([6]) one can find the ratio of reaction 
rates 

Jh 
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(7) 



In the case of Jh I Jnei.bg <C 1 the effect of neutral hydrogen 
on the Hell — > Hel recombination is negligible. In the case 

1 Following Seager et al. (2000) we consider in this section only 
2 1 p — » l 1 ^ transitions, although similar consideration can be pro- 
vided for Hel 2 3 p — > l 3 s transitions. 
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Figure 1. Hel level structure used in this paper 

of Jh / Jnei.bg 3> 1 effect of neutral hydrogen dominates and 
recombination speeds up essentially. 

Dependence of Jh I Jnei.bg on 2 is shown in fig. 
2. Exponential increase of Jh / Jnei.bg (with decrease of z) 
is concerned with exponential increase of neutral hydrogen 
concentration Nhi which is satisfied the Saha ionization 
equation for HI in the epochs z — 1500 — 3000. 

The ratio Jh / 'Jnei , bg is equal to I at epoch z ~ 1800 
(fig. 2) in contradiction with results by Seager et al. (2000) 
where the ratio Jh / JHei.bg is much less than unity in the 
period z = 1500 - 3000. 

Note that in paper by Boschan and Biltzinger (1998) 
other estimation of neutral hydrogen effect on helium recom- 
bination kinetics was proposed. In that work disappearance 
of Hel resonant photons was characterized by optical depth 
of hydrogen ionization by Hel 2 1 p — » I x s resonant photons 



r ~ ca H {v bg )NHi&.t 



(8) 



where At is the duration of period when considered photon 
is the Hel resonant photon. The value At is defined by de- 
viation between photon frequency and central frequency of 
Hel resonant transition because only Hel resonant photons 
keep the populations of Hel excited states. Using v — —Hv 
we find 



-^OH[Vb g )— — Nhi 

tl Vbg 

Thus one can obtain that 

Jh/ Jnel.bg 
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and we can see that both of estimation (by reaction rate 
ratio and by optical depth) give the same results. The value 
r is equal to 1 at z ~ 1800 according to (|I0l) (see fig. 2). 
This result is similar to the result of paper by Boschan and 
Biltzinger (1998). 

The simple estimation presented above shows that neu- 
tral hydrogen effect on helium recombination kinetics is es- 
sential. The fact made us consider this effect in detail and 
derive helium kinetic equation taking into account effect of 
neutral hydrogen. 
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3 MAIN EQUATION 

The behaviour of Hell (i.e. He + ) fraction is described by the following differential equation 



x = -C par \ a P arN e x - —j3 par exp I ~ j (1 - x ) j - Cor ya or N e x - -^-Por exp I --j^T£ J (1 - x )j (11) 

where x = Nue.il /N B e is the fraction of Hell ions relative to the total number of helium atoms and ions, C par is the factor 
by which the ordinary recombination rate is inhibited by the presence of Hel 2 1 p — > l 1 s resonance-line radiation, a par is the 
total Hell— >HeI recombination coefficient to the excited para-states of Hel, subscript a denotes state 2 1 s of Hel atom, N e is 
the free electron concentration, g a = 1 is the statistical weight of 2 1 s state of Hel, g g = 1 is the statistical weight of l 1 s state 
of Hel, [3 par is the total Hel^Hell ionization coefficient from the excited para-states of Hel, E ag is the 2 x s — *■ l 1 s transition 
energy, C or is the factor by which the ordinary recombination rate is inhibited by the presence of Hel 2 3 p — * Is resonance- line 
radiation, a or is the total Hell— >HeI recombination coefficient to the excited ortho-states of Hel, subscript a' denotes state 
2 3 s of Hel atom g a i = 3 is the statistical weight of 2 3 s state of Hel, (3 or is the total Hel— >HeII ionization coefficient from the 
excited ortho-states of Hel, E a t g is the 2 3 s — » l 1 s transition energy. 

The para- and ortho- recombination and ionization coefficients are related by the following formulae 

d V ar = ^a par g e (T) exp (-^] , Par = ^-a or g e (T) exp (12) 
g a \ kbI ) g a > V fc Bi / 

where subscript c denotes continuum state of Hel atom, g c = 4 is the statistical weight of continuum state of (He + +e~), 
g e (T) is the partition function of free electrons, E ca is the c — ► 2 1 s transition energy, E ca i is the c — > 2 3 s transition energy. 
The inhibition factor C par is given by the following formula 

c = (9b/ 'gg) {K + A r bg ) exp j~E ba /k B T) + A ag 

P P ar + {g b /ga){A^ + Al^e^ v {-E ba /k B T) + A ag 

where A bg is the effective coefficient [s _1 ] of 2 1 p «-» l 1 s transitions due to effect of neutral hydrogen on Hel 2 1 p — > l 1 s 
resonant radiation (superscript H denotes hydrogen), E ba is the 2 1 p — > 2 1 s transition energy, A ag is the coefficient of two- 
photon 2 1 s — *■ l 1 s spontaneous decay. 

The inhibition factor C or is given by the following formula 

(Wffa') (A$ g + Al, g ) exp (-E b , a ,/k B T) 

Cor — 7 r (14J 

Por + (gv/g a >) {A$ g + Al, g j exp (~E b , a ,/k B T) 

where subscript b' denotes state 2 3 p of Hel, g b > = 9 is the statistical weight of 2 3 p state of Hel, Ay is the effective coefficient 
of 2 3 p <-> f 1 s transitions due to effect of neutral hydrogen on Hel 2 3 p — > Is resonant radiation, A r b , g is the effective coefficient 
of 2 3 p «-> f 1 s transitions due to escape of photon from the 2 3 p — > l x s line profile because of cosmological redshift, E b > a i is the 
2 3 p — > 2 3 s transition energy. 

The coefficient A r bg is given by the formula ((5]). The coefficient Ay g is given by the following formulae 

il- a — l/i I \\ 9b' Ay g (? NheI ,g _v 

Ayg = 4{.' a V fl I 1 - ex P (- r 6' 9 )) » TVfl - (15) 

b'g 

where Ay g = g2 3 Pi^42 3 Pi^i 1 s /fli>' ' s the coefficient of 2 3 p — > l^s spontaneous transition, g2 3 Pi = 3 is the statistical weight 
of 2 3 Pi state. 

The coefficient Af g (subscript / denotes 2 1 p (b) or 2 3 p (b') state) can be obtained by consideration of kinetic equation 
for Hel / — > g transitions and transfer equation for Hel / — > g resonance radiation in the homogeneous expanding Universe in 
the presence of neutral hydrogen. Detailed calculation (paper in preparation) shows that exact expressions for the coefficients 
A bg and Ay can be approximated by the following formulae: 



,tf _ A fg = (9f/gg)A fg N He ic 2 

1 + P/7% ' V^th {u /g ) 8nu%Au Difg N HI 



A fg= , — jr f > 1f9= 7^Z7^. ^ ^■2 A ,._ . at... ( 16 ) 



Parameters pf and qf are the following: p b = 0.36, q b = 0.97, p b < = 0.66, q b < = 0.9. 

The equation (|lf | has been solved numerically. Results of calculation are presented in fig. 3. 

The equation combined with equations for description of fractions of HII (Zel'dovich et al., 1968; Peebles, 1968) 
and Helll (Matsuda et al., 1969; Seager et al., 1999) allows us to calculate ionization history of the Universe in the epochs 
2 = 10 2 - 10 4 . 
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Figure 2. The dependence of the ratio Jh/ J Hel.bg on redshift 
z is shown. The dotted curve shows the unity level. 



4 RADIATION CONCERNED WITH 
HEII^HEI RECOMBINATION 

The most of recombination photons corresponding to tran- 
sitions to ground state of helium atoms is absorbed by the 
neutral hydrogen atoms and converts to the Lya-photons of 
HI. 

Number of HI Lya-photons per mode concerned with 
Hell— >HeI recombination is given by the formula 



Arj(v a , z b ) = 



c 3 7V ge (z b ) \x(z b )\ 
8tt^ H{z b ) 



(17) 



where Arj is the number of photons per mode depending on 
frequency v and epoch z, v a is the HI Lya frequency, z b is 
the epochs of photon birth. Using relation 



Ar]{v, z) — Arj ( u- ^ 



\ 1 + z 



(18) 



we can find the spectrum of radiation concerned with 
Hell— >HeI recombination at the present epoch. 



5 RESULTS 

The main results of this paper are the dependencies of Hell 
fraction on redshift (fig. 3) and free electron fraction on 
redshift (fig. 4). The calculation shows that during epochs 
z — 1600 — 2000 kinetics of Hell— >HeI recombination (this 
work, black curve) changes from strongly non-equilibrium 
(Seager et al., 1999, blue curve) to quasi-equilibrium (ac- 
cording to the Saha formula for Hel, red curve). This change 
is due to absorption of Hel resonant photons by neutral hy- 
drogen whose concentration increases exponentially with de- 
crease of temperature. At the epoch z — 1600 the fraction of 
Hell relative to the total number of helium ions and atoms is 
less than 10" 9 . To compare our results with those of previous 
papers we have calculated Hell^Hel recombination kinetics 
at the following parameters: Af g — 0, A 2 3p 1 ^ 1 i g = 233 s _1 
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Figure 3. The dependence of Hell fraction x on redshift z for var- 
ious model of Hell— >HeI recombination is presented: blue curve 
corresponds to result of Seager et al., 1999, violet curve corre- 
sponds to taking into account recombination through ortho-states 



(Dubrovich and Grachev, 2005) and j4 2 3 p 1 ^i 1 



So 



233 s _1 



(Lin et al., 1977), magenta curve corresponds to taking into ac- 
count recombination through ortho-states (Wong and Scott, 2006) 
and A 2 3p 1 _ >1 i So = 177.58 s -1 (Lach and Pachucki, 2001) [vio- 
let curve practically overlaps magenta curve], green curve cor- 
responds to taking into account effect of neutral hydrogen on 
Hell— »HeI recombination, black curve corresponds to taking into 
account both of the effects, red curve corresponds to recombina- 
tion according to the Saha formula for Hel. Top panel corresponds 
to linear scale, bottom one does logarithmic scale. 



(Lin et al., 1977) corresponding to paper by Dubrovich and 
Grachev (2005) - violet curve, and ^4 2 3 p 1 -»i 1 s ~ 177.58 s~ 
(Lach and Pachucki, 2001) corresponding to paper by Wong 
and Scott (2006) - magenta curve. These results show very 
good agreement with results by Dubrovich and Grachev 
(2005) and Wong and Scott (2006). 

In the fig. 4 one can see that recombination of helium 
ends (x ^ 10~ 6 ) before recombination of hydrogen begins 
(xhi S3 0.99) in contradiction of results by Seager et al. 
(1999). This change of primordial plasma recombination ki- 
netics should led to changes of calculated CMBR anisotropy 
at the level of up to 5% for multipoles I ~ 1500 — 3000 
(Seager et al., 2000). 

Results of calculation of radiation intensity are pre- 
sented in the fig. 5: 1) Planck spectrum radiation at the 
temperature To = 2.726 K (red curve); 2) HI Lya radia- 
tion concerned with hydrogen recombination (blue curve, 
Grachev and Dubrovich, 1991; Rubino-Martin et al., 2006; 
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Figure 4. The dependence of free electron fraction on redshift z 
for various model of Hell— >HeI recombination is presented: blue 
curve corresponds to result of Seager et al., 1999, black curve 
corresponds this work, red curve corresponds to recombination 
according to the Saha formula. 

Kholupenko and Ivanchik, 2006; and references therein); 
3) HI Lya radiation concerned with Hell— >HeI recombi- 
nation (magenta curve). Wien-tail CMBR intensity dis- 
tortion concerned with Hell— >HeI recombination has the 
maximum value about 9.3 • 10 -25 erg-cm _2 ster _1 s _1 Hz _1 
being at frequency 1350 GHz. This frequency is about 2 
times smaller than frequency calculated by Wong et al. 
(2006). This difference is concerned with re-emission of 
helium recombination photons by neutral hydrogen atoms 
in Lya-line of HI while in the paper by Wong et al. 
(2006) helium recombination photon emission was consid- 
ered at transitions 2 1 p — > l x s and 2 1 s — > l 1 s of Hel. 



Table 1. Parameters of 
model 


the standard 


cosmological 


Value description 


Symbol 


Value 


total matter 


fitot 


1 


non-relativistic matter 




0.27 


baryonic matter 


fib 


0.04 


relativistic matter 




~ io- 4 


vacuum-like energy 


fi A 


0.73 


Hubble constant 


H 


70 km/s/Mpc 


radiation temperature 


To 


2.726 K 


helium mass fraction 


Y 


0.24 
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v, GHz 

Figure 5. The dependence of radiation intensity on frequency 
at present epoch: red curve corresponds to Planck radiation at 
temperature To = 2.726 K, blue curve corresponds to HI Lyo ra- 
diation concerned with recombination of hydrogen, magenta curve 
corresponds to HI Lyo radiation concerned with recombination 
of helium (HeII->HeI) 
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